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Abstract. - Recent measurements of the temperature fluctuations of the Cosmic Microwave 
Background (CMB) radiation from the WMAP satellite provide indication of a non-Gaussian 
behavior. Although the observed feature is small, it is detectable and analyzable. Indeed, 
the temperature distribution P CMB (AT) of these data can be quite well fitted by the anoma- 
lous probability distribution emerging within nonextensive statistical mechanics, based on the 
entropy S q = k{l - J dx [P{x)] q }/(q - 1) (Si = -k J dx P(x)ln[P(x)]). For the CMB fre- 
quencies analysed, v = 40.7, 60.8, and 93.5 GHz, P CMB (AT) is well described by P 9 (AT) cx 
1/[1 + (q- l)B(>)(AT) 2 ] 1/(9-1) , w ith q = 1.04 ± 0.01, the strongest non-Gaussian contribution 
coming from the South-East sector of the celestial sphere. Moreover, Monte Carlo simulations 
exclude, at the 99% confidence level, Pi (AT) cx e - s M( AT ) to fit the three -year data. 



. ^ Statistical properties of the CMB temperature 
fluctuations. — Since the early 90's a set of pre- 
5h else observations showed that the angular distribution 



of the Cosmic Microwave Background (CMB) radiation 
presents tiny fluctuations, of the order of 10 -5 , around 
the equilibrium temperature T . These CMB tempera- 
ture fluctuations, denoted AT(#, c/>), are functions of the 
direction of observation in the celestial sphere (8, </>) and 
can be expressed by using a spherical harmonic expansion 
AT(0,<p)/T Q = J2e m ae m Y lm (6,(j>). Standard cosmologi- 
cal theory predicts that the complex ag rn modes are statis- 
tically isotropic and Gaussian random fields of zero mean 
((dim) = 0) and non-zero variance ((a^, , a tm ) ^ 0). The 
fact that AT is a real quantity yields consistency rela- 
tionships between the real and imaginary parts of ai m 
and at.- m . Gaussian randomness means that both real 
and imaginary parts of the a£ m are each an independent 
random variable obeying a Gaussian distribution of zero 
mean. Statistical isotropy means that the variance, which 
in principle could be different for each £ and m, depends 
only on I, i.e. (a}, m ,a tm ) = Ct5yt5 m > m . Thus, in the 
standard inflationary big-bang model, one expects that 
the CMB temperature fluctuations AT should be appro- 



priately described by a Gaussian distribution. 

The recent release of the three-year high angular reso- 
lution and low instrument noise CMB measurements from 
the Wilkinson Microwave Anisotropy Probe (WMAP) [1] 
offered the opportunity to analyse, with unprecedent pre- 
cision, their temperature distribution in order to test 
the Gaussian hypothesis. These WMAP maps [2], with 
parameter resolution N S id e = 512, divided the celestial 
sphere into 3, 145, 728 equal sized areas, termed pixels, 
each one with an assigned highly precise AT value. Be- 
cause measurements in the region around the Galactic 
Plane are strongly contaminated with foregrounds, they 
shall be eliminated from the data set in analysis by apply- 
ing a cut-sky mask. The remaining ~ 2.5 million pixels 
guarantee a fine adjustment of the temperature distribu- 
tion histogram (also termed one-point distribution func- 
tion) to a smooth curve, thus minimizing the uncertainty 
of any (Gaussian or non-Gaussian) fitting. 

Through the Minkowski functional's analysis, the 
WMAP team placed tight new limits on non-Gaussian 
temperature fluctuations in the CMB, thus concluding 
that these data are consistent with Gaussianity [3, 5]. 
However, several kinds of non- Gaussian features have been 
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reported in the first-year and three-year WMAP data [6] . 
Clearly, the detailed study of the Gaussian random hy- 
pothesis must take into account the possibility that devia- 
tions from Gaussianity may have non-cosmological origins, 
such as unsubtracted foreground contamination, instru- 
mental noise, and/or systematic effects [7]. Nonetheless, it 
is also possible that deviations from Gaussianity, although 
small, may be intrinsic to the CMB radiation [8] or have 
some other cosmological origin such as cosmic strings [9] . 
The related issue regarding the violation of the large-scale 
statistical isotropy, termed the North-South asymmetry, 
has also been extensively investigated in a variety of anal- 
yses, using different mathematical-statistical tools, with 
the first-year WMAP data [10] and more recently also 
with the three-year data [11]. 

Given the potential importance of these reported non- 
Gaussian features for our understanding of the physics of 
the early universe, all these analyses must be considered as 
non-exhaustive and further studies have to be performed. 
Although these non-Gaussian findings [8] - obtained with 
the first-year WMAP data - could be interpreted as un- 
removed foregrounds and/or instrumental noise, a com- 
parative analysis with three-year data, where the signal 
to noise ratio increases due to a greater number of obser- 
vations per pixel, could confirm or discard some of these 
possible explanations, or simply put better limits of their 
effect on the CMB data. This converts the analyses of the 
statistical properties of the WMAP three- year data into a 
necessary and important task. 

Nonextensive statistical mechanics appears to be suit- 
able to describe statistical properties of long-range cor- 
related systems [12], like the primordial CMB radiation 
filling the universe. For such systems, correlations are 
expected to break down the Gaussian, i.e., essentially un- 
correlated, properties [13]. In that situation, the q-index 
can be seen as a scalar measure of the correlation range 
present in the system. For example, for an attractive two- 
body potential decaying as l/r a , a > 0, in a d-dimensional 
classical system, we may expect q to increase above unity 
when a/d decreases below unity [14]. 

In what follows we perform a brief overview of the non- 
Gaussian distributions emerging from nonextensive statis- 
tical mechanics. 

Basics of nonextensive statistical mechanics. 

The probability distribution function P q in the frame 
of the nonextensive statistical mechanics results [8] from 
the optimization of the g-entropy defined [15] (see [16] 
for various applications, and [17] for a new experimental 
confirmation) as follows (q £ 3?) 



S 1 = S BG = -kJ dxP(x) In[P(x)] , 

satisfying suitable constraints (BG stands for Boltzmann- 
Gibbs). Given two independent systems IU,Eb ( m the 



sense that P^aUVb = pS A pS B ^ tlle pseu doadditivity 
property says that Sf AUllB = Sf A +Sf B + (l-q)Sf A Sf B . 
In other words, (1 — q) is a measure of the nonadditivity of 
the total system XUUSb. Clearly, the standard, extensive, 
statistical mechanics is recovered in the q — > 1 limit. 
From Eq. (1) one deduces [8] (see also [18] regarding the 
connection with the g-Central Limit Theorem) 

P q (AT)=A q e- B ^\ (2) 

where in general B q depends on the q-index, and A q is 
the normalization constant obtained in such a way that 
J dx P q (x) = 1. The ^-exponential is defined by 

e\ = [1 + (1 - q)z] 1 '^ , for [1 + (1 - > 0, (3) 

and e q = otherwise. In the limit q — > 1, we recover the 
Gaussian distribution 

Pq(AT) -» Pi = pGauss( AT ) =Ae -B(u)(ATf ^ ^ 

where A = l/(er„\/27r), B{v) = 1/(2ct 2 ), and er 2 , is the 
variance of the Gaussian distribution. 

We shall use the nonextensive probability distribution 
P q , defined in Eqs. (2) and (3), to study the CMB tem- 
perature distribution p CMB from the WMAP three-year 
data. Notice that by plotting the temperature distribu- 
tion histogram in the form of the Number of Pixels ver- 
sus (AT/cr„) 2 any Gaussian distribution fits a straight 
line, thus recognizing departures from Gaussianity (i.e., 
from linearity) becomes easy. For a clear comparison be- 
tween the Gaussian (linear) and the nonextensive distri- 
butions we shall plot both distributions P q and pGauss 
having equal value at the initial point AT = (that is, 
A q = A = 1). Moreover, we observe that the best-fitting 
tangent line at that point has equal slope for both distri- 
butions, which implies that B q = B(v) (this procedure is 
consistent with the fact that deq/dz\ z= a = 1, Vq). There- 
fore, the best-fitting of the CMB data through the nonex- 
tensive distributions is obtained by adjusting only the q- 
index. 

WMAP data analyses. - We analyse the WMAP 
three-year data by exploring their temperature distribu- 
tion in terms of Gaussian and nonextensive probability 
distributions. Since any non-Gaussian signature is re- 
vealed as non-linearities when the temperature distribu- 
tions are plotted in the form of the Number of Pixels ver- 
sus (AT/ov) 2 , we present our results using this type of 
plots. We use the x 2 /degree of freedom (dof) estimator 
test to fit the data with the best Gaussian and nonexten- 
sive distributions functions. 

We investigate these data corresponding to the Q, 
V, and M^-bands (central frequencies at 40.7, 60.8, and 
93.5 GHz, respectively) by analysing the eight foreground 
reduced differencing assemblies (DAs) available at [2], 
namely Ql, Q2, VI, V2, Wl, W2, W3, and WA. These 
DA CMB maps were corrected by the WMAP team for 
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the Galactic foregrounds (free-free, synchrotron, and dust 
emission) using the 3-band, 5-parameter template fitting 
method described in [1]. However, the foreground removal 
is only applicable to regions outside the Kp2 mask (which 
cuts ~ 15% of the sky data). The more severe cut, rep- 
resented by the KpO mask (which excludes ~ 23% of the 
sky data), is recommended [1] for the statistical analysis 
of the CMB maps. Thus, for our analyses we have used 
the KpO mask. For the sake of completeness, we have 
also examined the foreground-reduced coadded maps Q (= 
Q1®Q2), V(= VI® V2), and W (= Wl® W2® W3® W4), 
as well as the coadded V® IF map, formed by the six DAs 
VI® V2® Wl® W2® W3® WA, and the coadded Q® V® W 
map, formed by the eight DAs (the ® symbol means the 
coadded sum of CMB maps as defined in [19]). 

Finally, we also investigated the temperature distribu- 
tion in the four quadrants of the celestial sphere (NE, 
NW, SW, and SE) in order to test the claimed large-scale 
anomaly [10, 11] in the CMB temperature fluctuations 
data. 

As pointed out in [9] , the signal measured at any pixel in 
the microwave sky is made of several components Tpi xe i — 
T foregrounds+ T noise+ T CMB, corresponding to foreground 
contamination signals, the noise from the instruments, and 
the CMB temperature fluctuations, respectively. Fore- 
ground contributions are expected around the Galactic 
Plane and from point sources. Given the fact that the 
instrument noise in the map is inhomogeneous (higher in 
the Ecliptic Plane than near the poles), the contribution 
from T no i se is important. For the subtraction of the fore- 
grounds from CMB maps we used, as mentioned above, 
the most severe cut-sky, i.e., the KpO mask (for compari- 
son, we also made the analyses with the Kp2 mask). The 
masks are provided by the WMAP team [1] and allow for 
the selective exclusion of portions of the sky basically due 
to contaminating radiative processes from our galaxy, in- 
cluding a 0.6 degree radius exclusion area around known 
point sources in the celestial sphere. The application of 
a mask to a CMB map simply means that, for a given 
pixel, a mask value of zero implies that it is excluded; a 
value of one means that it is accepted. Moreover, since the 
Galactic foregrounds contribute as positive temperature 
contaminations, we perform the analyses of the one-point 
distributions considering only the negative temperature 
fluctuations data. In this way, the first term in the above 
equation is dropped and one is left with CMB signal plus 
the signal from the instrument noise (the signal noise is 
Gaussian per observation [9,19]). What converts the in- 
strument noise in a non-Gaussian signal is the non-uniform 
way in which the CMB sky was observed, because some 
regions were largely more observed than others, making 
such a noise to be different from pixel to pixel. The way 
out to deal with this problem is to characterize its shape 
and intensity through Monte Carlo realizations of similar 
sky data observations. 

Consider the remaining effects by defining [9] the vari- 



ance of the total signal Tp ixe j as 



a » = CT pixel 



^0, 2 

m CMB 



(5) 



where n, is the number of observations for the ith pixel, 
<7q is the noise variance per observation of a given instru- 
ment [1], and Oq MB is the variance of the CMB temper- 
ature fluctuations. The mean contribution of the instru- 
mental noise can be estimated by considering the effect of 
the different number of observations for each pixel. This 
can be done with the effective variance due to the non- 
stationary instrument noise [9] 
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N 



(6) 



where N is the total number of pixels considered in the 
analysis of the map. The values of XXV n i) f° r each 
WMAP DA map are given in table 1 (KpO mask case) 
The total number of pixels N with CMB signal corrected 
for Galactic contaminations results from the application 
of the cut-sky mask, that is N = N Kp0 = 2,407, 737 for 
the KpO mask, or N = N Kp2 = 2,664,479 for the Kp2 
mask. Thus, the effective noise variances cr ^ o i S0 (^) and 
the variances of, leads to the CMB variance 



CMB 



>)• 



(7) 



Although both ct 2 and c^oise^) depend on the map un- 
der analysis, their difference is independent of the map. 
In other words, if the treatment of instrumental noise is 
correct, the CMB variance cr^MB should be the same for 
the maps under investigation (Q, V, and W). 

The present best-fit data analyses were obtained ac- 
cording to the x 2 /degree of freedom (dof) estimator test. 
Thus, the typical x 2 /200 values for the Gaussian fits are 
2 x 10~ 2 , and the x 2 /200 values for the nonextensive tem- 
perature distribution arc 2 x 10~ 3 . Analyses performed 
using 300 dof and 400 dof, instead of 200 dof, result in sim- 
ilar best-fitting values. In addition, we have not detected 
any significative difference using the Kp2 mask instead of 
the KpO mask. 

Our analyses of the eight DAs is as follows. Once the 
variances ct 2 have been determined through the x 2 best- 
fit Gaussian temperature distribution (blue lines in fig. 1), 
we use the effective noise variance &noise'> §i ven m Eq- (6), 
to calculate the CMB variances o-q MB ( see table 1). This 
approach validates, although a posteriori, the use of the 
effective noise variance as representing the mean contri- 
bution of the instrumental noise and our previous results. 
The CMB variances we obtain for the Q, V, and W coad- 
ded maps are cr^ MB = 5.98, 6.56, 6.77 x 10 3 ^K 2 , respec- 
tively. 

Summarizing, the temperature distribution analyses of 
the coadded Q, V, and W CMB maps, as well as the eight 
individual DAs maps shown in fig. 1, exibit that the dis- 
tribution of the CMB temperature fluctuations does not 
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Fig. 1: Top: Fits to the (positive and negative) WMAP CMB temperature fluctuations data, corresponding to the Q, V, and 
coadded maps (after the KpO cut-sky), in the number of pixels versus AT/cr„ plots. We show the \ 2 best-fits: Gaussian 
distribution (blue curve) with <tq = 104 fiK, ay — 118 fiK, and cr\y = 131 fiK, respectively, and each nonextensive distribution 
(red curve) P q with q = 1.04. Bottom: Similar analysis, but now with the eight DA WMAP maps (Q1,...,W4) after applying 
the KpO mask. We plotted the number of pixels versus (AT/a v ) 2 to enhance the non-Gaussian behavior. To avoid possible 
unremoved Galactic foregrounds, we consider only the negative temperature fluctuations. Again, we show the x best-fits: 
Gaussian distribution (blue curve) and nonextensive distribution (red curve) P q , now with q — 1.04 ± 0.01. 



obey a Gaussian distribution. The analysis of the other 
coadded CMB maps, mentioned before, produces a result 
that is fully consistent with this one. 

Another way to test Gaussianity in a given data set is 
the calculation of the factor involving the 4th- and 2nd- 
momcnta of the distribution, R = (AT 4 ) / (3AT 2 ) 2 , where 
(x n ) = J dxP(x)x n (R=l corresponds to the Gaussian 
distribution case). All the CMB maps analysed results in 
R € [1-02, 1.06], which validates the possible description, 
within the four available ordinate decades, of the WMAP 
data using the P q distributions emerging from nonexten- 
sive statistical mechanics. 

We also studied the coadded W map, after applying 
the KpO mask, in the four quadrants: North-East (NE), 
North- West (NW), South- West (SW), and South-East 
(SE) of the celestial sphere. It turns out that the tem- 
perature distribution of the CMB data is sensibly differ- 
ent in one of these four sky patches. Indeed, we found 

q NE - L015 ' q NW ~ L01 ' q SW ~ L02 ' and %E ~ 1 - 05 ' 

with R = 1.02, 1.01, 1.025, and 1.04, respectively. 

The significance of the g-index values found are based 



on the robustness of our findings. In fact, all our results 
are the same independently from the frequency-map or 
the coadded-map analysed. Also, they are equal for the 
two cut-sky masks applied, and for the one-year [8] and 
three-year maps. 

Monte Carlo analysis. — Now we present the results 
of our analyses considering the WMAP Monte Carlo (MC) 
simulations, which contain the instrument noise having a 
spatially varying noise structure, like real data. 

As a matter of fact, the WMAP detectors noises are 
inhomogeneous. This is due to the fact that the re- 
gions around the ecliptic poles (NW and SE quadrants) 
have large number-of-observations (A^g), and therefore 
low statistical noise, while the region corresponding to 
the ecliptic plane (NE and SW quadrants) were less ob- 
served (lower iV 0DS values), and have large statistical noise. 
The spatially varying noise structure for each DA CMB 
map, with simulated instrument noise that include all 
known radiometric effects, were properly realized by the 
WMAP team in a set of MC sky realizations for each 
DA (available at http:\\lambda.gsfc.nasa.gov\product\ 
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map\drl\noise_sims_imaps.cfm). These MC data con- 
stitute an important tool for instrument noise evaluation 
corresponding to each DA, and also for the MC-coadded 
maps, obtainable for each one of the three bands (Q, V, 
and W) by combining MCs from different DAs. For in- 
stance, to obtain a MC-coadded CMB map in the Q band, 
one chooses QV from the set of MCs for the DA-(J1, and 
combine (in a coadded way, according to [19]) with Q2- 7 
from the set of MCs for the DA- Q2. 

This time, the noise instrument of a given map is treated 
by computing the variance-normalized temperature for 
each pixel (a procedure described in [4]: see, in partic- 
ular, Eq. (17)). We observe that the instrument noise 
shows non-Gaussian behavior whose strength and signa- 
ture are different to what is observed in our analyses of 
the 3-year WMAP data. In average, we obtained \ 2 best- 
fits with distributions P q with q = 1.1 ± 0.02, and after 
the variance-normalization, the g-index maintain its value, 
although the variance is clearly modified. 
Moreover, we noticed that the non-Gaussian features of 
the instrument noise are sensible to N ^ s , in the sense 
that the g-index of the P q distributions changes signifi- 
cantly when analysing two regions in the sky, namely, one 
that was observed (at least) twice the other. This inter- 
esting (and expected) feature is absent in the three-year 
WMAP data when compared with the (less observed) one- 
year WMAP data [8], in fact we obtained essentially the 
same g-index for both data sets. Finally, we have also to 
mention that such strong non-Gaussianities, present in all 
the coadded-MCs, exhibit a high symmetry between the 
SE and NW quadrants: <Z SE /<Z NW = 1.0 ± 10~ 3 . Again, 
this feature is fully absent in WMAP data, where we have 

«se/«nw ~ lM - 

Comments and Conclusions. — Some comments 
are in due regarding the use of DA maps and the applica- 
tion of two WMAP masks, KpO and Kp2, in our analysis. 
First, we analysed frequency by frequency (Q, V, and W 
bands) because the CMB Galactic foregrounds (free-free, 
synchrotron, and dust emission) are frequency dependent, 
which means that, if they are still present in the WMAP 
data, they should manifest differently at each frequency. 
Moreover, as recommended by the WMAP team, we have 
analysed the foreground reduced CMB maps [2], which 
were produced using the Foreground Template Model to 
subtract the synchrotron, free-free, and dust emission from 
the measured sky maps [1]. Our analyses of the one- 
point distributions of all these WMAP CMB maps do 
not show significant differences in the value of the g-index 
for the data investigated. Second, another way to exam- 
ine the possible presence of unremoved foregrounds is by 
analysing the CMB data using different cut-sky masks. 
We performed such analyses using both the KpO and Kp2 
masks and, according to the \ 2 best-fit estimator, we do 
not found significant difference in the values of the g-index. 
Third, known Galactic foregrounds contribute with posi- 
tive temperature fluctuations contaminations. Hence, we 



consider only the negative temperature fluctuations of the 
CMB maps to perform our data analyses (see the bottom 
plots of fig. 1). These procedures, and their correspond- 
ing results, show that possible unremoved Galactic fore- 
grounds could not be responsible for the detected non- 
Gaussian signatures. 

Now consider the instrument noise, which is always con- 
volved with the CMB signal. In WMAP maps the instru- 
ment noise is inhomogeneous, higher in the Ecliptic Plane 
(roughly NE and SW quadrants) than near the Ecliptic 
Poles (NW and SE quadrants), and their contribution to 
-^pixel ma y cause a non-Gaussian behavior, as one can ob- 
serve in MC sky maps simulated by the WMAP team. 
However, two interesting facts make the non-Gaussian 
noise signals qualitative and quantitatively different from 
what is observed in WMAP data. We emphasize that, 
in order to deal with inhomogeneous instrument noise in 
our CMB analyses, we have taked into account the fol- 
lowing two procedures: the coadded sum of CMB DAs 
maps, as described in [19], and the variance-normalized 
temperature, as described in [4]. First, the instrument 
noise is extremely sensible to N Q ^ S , as expected, and as 
MC simulations confirm. It is well known that sky pixels 
in three-year WMAP data has been observed roughly 3 
times more than one-year data, nevertheless the P q dis- 
tributions are the best-fit one-point dsitributions, with 
q = 1.04 ± 0.01, for both data sets (see [8]). Second, 
the non-Gaussian signals from inhomogeneous instrument 
noise exibit a high symmetry between the NW and SE 
quadrants that has been quantifiable using the MC maps: 
g gE /g NW = 1.0 ± 10~ 3 , nonetheless this feature is fully 
absent in WMAP data, where we have Q se /q nw — 1-04. 
These two facts strongly indicate that any subtle effect 
coming from inhomogeneous instrument noise contamina- 
tion could not explain the presence of non-Gaussian sig- 
natures in the analysed WMAP data. 

Summarizing, we have shown that a Gaussian distri- 
bution is excluded, at the 99% CL, to properly represent 
the overall CMB temperature fluctuations measured by 
WMAP. Through the temperature distribution analyses 
we notice that the three-year WMAP data are suitably 
explained, along about four decades, by a P q distribution 
with q = 1.04 ±0.01. Moreover, our results also evidence 
that the strongest non-Gaussian contribution comes from 
the South-East quadrant, where g gE ~ 1.05. Therefore, 
the hypothesis of a nonextensive nature for the CMB tem- 
perature fluctuations is quite plausible. 
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Table 1: All the variances are in units: xl0 3 ^K 2 . Analyses of the 8 DA CMB maps, A?KpO = 2,407,737 pixels. DE means 
'direct evaluation' from the corresponding CMB map. According to the x 2 best-fit estimator, the temperature distributions of 
these CMB data are described by P q distributions with variance ol and g-index: q — 1.04 ± 0.01. 
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